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Massless Dirac electrons in graphene fill Landau levels with energies scaled as square roots of their 
numbers. Coulomb interaction between electrons leads to mixing of different Landau levels. The 
relative strength of this interaction depends only on dielectric susceptibility of surrounding medium 
and can be large in suspended graphene. We consider influence of Landau level mixing on the prop- 
erties of magnetoexcitons and magnetoplasmons — elementary electron-hole excitations in graphene 
in quantizing magnetic field. We show that, at small enough background dielectric screening, the 
mixing leads to very essential change of magnetoexciton and magnetoplasmon dispersion laws in 
comparison with the lowest Landau level approximation. 

PACS numbers: 73.22.Pr, 71.35.Ji, 73.43.Mp, 71.70.Gm 



INTRODUCTION 

Two-dimensional systems in strong magnetic field 
arc studied intensively since the discovery of integer 
and fractional quantum Hall effects [1-3]. For a long 
time, such systems were represented by gallium arsenide 
heterostructures with two-dimensional electron motion 
within each subband [4]. 

New and very interesting realization of two- 
dimensional electron system appeared when graphene, 
a monoatomic layer of carbon, was successfully isolated 
[5, 6]. The most spectacular property of graphene is the 
fact that its electrons behave as massless chiral particles, 
obeying Dirac equation. Intensive experimental and 
theoretical studies of this material over several recent 
years yielded a plethora of interesting results [7-9]. In 
particular, graphene demonstrates unusual half-integer 
quantum Hall effect [6], which can be observed even at 
room temperature [10]. 

In external perpendicular magnetic field, the mo- 
tion of electrons along cyclotron orbits acquires zero- 
dimensional character and, as a result, electrons fill dis- 
crete Landau levels [11]. In semiconductor quantum 
wells, Landau levels are equidistant and separation be- 
tween them is determined by the cyclotron frequency 
ojc = eH/mc. In graphene, due to massless nature of elec- 
trons, "ultra-relativistic" Landau levels appear, which 
are non-equidistant and located symmetrically astride 
the Dirac point [12, 13]. Energies of these levels are 
En = sign{n)y/2\n\vF /Ih, where n ~ 0,±1,±2,..., 
vp ~ 10^ m/s is the Fermi velocity of electrons and 
Ih = \/ c/eH is magnetic length, or radius of the cy- 
clotron orbit (here and below we assume ft = 1). 

In the case of integer filling, when several Landau levels 
are completely filled by electrons and all higher levels are 
empty, elementary excitations in the system are caused 
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by electron transitions from one of the filled Landau lev- 
els to one of the empty levels [14]. Such transitions can 
be observed in cyclotron resonance or Raman scatter- 
ing experiments as absorption peaks at ceratin energies. 
With neglect of Coulomb interaction, energy of the ex- 
cited electron-hole pair is just a distance between Landau 
levels of electron and hole. Coulomb interaction leads to 
mixing of transitions between different pairs of Landau 
levels, changing the resulting energies of elementary ex- 
citations. 

Characteristic energy of Coulomb interaction in mag- 
netic field is /eln, where e is a dielectric permittivity of 
surrounding medium. The relative strength of Coulomb 
interaction can be estimated as ratio of its characteristic 
value to a characteristic distance between Landau levels. 
For massive electrons in semiconductor quantum wells, 
this ratio is proportional to H~^/^, thus in asymptoti- 
cally strong magnetic field the Coulomb interaction be- 
comes a weak perturbation [15, !(>]. In this case, the low- 
est Landau level approximation, neglecting Landau level 
mixing, is often used. It was shown that Bose-condensate 
of noninteracting magnetoexcitons in the lowest Landau 
level as an exact ground state in semiconductor quantum 
well in strong magnetic field [17]. 

A different situation arises in graphene. The relative 
strength of Coulomb interaction in this system can be 
expressed as Vg = e'^/ev-p and does not depend on mag- 
netic field [in]. The only parameter which can infiuence 
it is the dielectric permittivity of surrounding medium 
£. At small enough e, mixing between different Landau 
levels can significantly change properties of elementary 
excitations in graphene. 

Coulomb interaction leads to appearance of two types 
of elementary excitations from the filled Landau levels. 
From summation of "ladder" diagrams we get magne- 
toexcitons, which can be imagined as bound states of 
electron and hole in magnetic field [14, 16, 19]. Prop- 
erties of magnetoexcitons in graphene were considered 
in several works, mainly in the lowest Landau level ap- 
proximation [20-24]. At e w 3. Landau level mixing was 
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shown to be weak in the works [20, 25]. 



I. MAGNETOEXCITONS 



Note that influence of Landau level mixing on proper- 
ties of an insulating ground state of neutral graphene was 
considered in by means of tight-binding Hartrce-Fock 
approximation. It was shown that Landau level mixing 
favors formation of insulating charge-density wave state 
instead of ferromagnetic and spin-density wave states in 
suspended graphene, i.e. at weak enough background di- 
electric screening. 

From the experimental point of view, the most inter- 
esting are magnetoexcitons with zero total momentum, 
which are only able to couple with electromagnetic radia- 
tion due to very small photon momentum. For usual non- 
relativistic electrons, magnetoexciton energy at zero mo- 
mentum is protected against corrections due to electron 
interactions by the Kohn theorem [27]. However, for elec- 
trons with linear dispersion in graphene the Kohn theo- 
rem is not applicable [21, 24, 28-32]. Thus, observable 
energies of excitonic spectral lines can be seriously rcnor- 
malized relatively to the bare values, calculated without 
taking into account Coulomb interaction. 

The other type of excitations can be derived using 
the random phase approximation, corresponding to sum- 
mation of "bubble" diagrams. These excitations, called 
magnetoplasmons, are analogue of plasmons and have 
been studied both in two-dimensional electron gas [14, 33] 
and graphene [18, 20, 21, 24, 34-39] (both with and with- 
out taking into account Landau level mixing). 

In the present article, we consider magnetoexcitons 
and magnetoplasmons with taking into account Landau 
level mixing and show how the properties of these ex- 
citations change in comparison with the lowest Landau 
level approximation. For magnetoexcitons, we take into 
account the mixing of asymptotically large number of 
Landau levels and find the limiting values of cyclotron 
resonance energies. 

For simplicity and in order to stress the role of vir- 
tual transitions between different pairs of electron and 
hole Landau levels (i.e. the role of two-particle pro- 
cesses), here we do not take into account renormalization 
of single-particle energies via exchange with the filled lev- 
els. This issue have been studied in several theoretical 
works [20, 21, 24, 30, 40]. Correction of Landau level 
energies can be treated as renormalization of the Fermi 
velocity, dependent on the ultraviolet cutoff for a number 
of the filled Landau levels taken into account in exchange 
processes. 

The rest of this article is organized as follows. In Sec- 
tion 2 we present a formalism for description of magne- 
toexcitons in graphene, which is applied in Section 3 to 
study influence of Coulomb interaction and Landau level 
mixing on their properties. In Section 4 we study mag- 
netoplasmons in graphene in the random phase approxi- 
mation and in Section 5 we formulate the conclusions. 



Electrons in graphene populate vicinities of two 
nonequivalcnt Dirac points in the Brillouin zone, or two 
valleys K and K'. We do not consider intervalley scat- 
tering and neglect valley splitting, thus it is sufficient to 
consider electrons in only one valley and treat existence 
of the other valley as additional twofold degeneracy of 
electron states. 

We consider magnetoexciton as an electron-hole pair, 
and we will denote all electron and hole variables by the 
indices 1 and 2 respectively. In the valley K, Hamilto- 
nian of free electrons in graphene in the basis {Ai,B2} 
of sublatticcs takes a form [7]: 



HI 



(0) 



(1) 



where pi± = {pix ± ipiy)/V2 are the cyclic components 
of electron momentum and wp ~ lO^m/s is the Fermi 
velocity of electrons. 

For external magnetic field H. parallel to the z axis, 
we take the symmetrical gauge, when A(r) = i[H x 
r]. Introducing the magnetic field as substitution of the 
momentum Pi — > Pi + (e/c)A(ri) in (1) (we treat the 
electron charge as — e), we get the Hamiltonian of the 
form: 



Hi 



V2 



H 



ai 
at 



(2) 



Here the operators ai ~ IrPi- — i?'i-/2ZH and = 
IhPi+ + iri+/2lH (where ri± = {xi ± iyi)/V2) obey 
bosonic commutation relation [ai,a^] = 1. 

Using this relation, by means of successive action of 
the raising operator we can construct Landau levels 
for electron [18] with energies 



and wave functions 



Sn(f>\n\-l,k(j) 
0|rx|fc(r) 



(3) 



(4) 



Here fc — 0, 1, 2, . . . is the index of guiding center, which 
enumerates electron states on the n-th Landau level {n = 
—00, . . . , -l-oo), having macroscopically large degeneracy 
= equal to a number of magnetic flux quanta 

penetrating the system of the area S. Eigenfunctions 
(pnk (r) of a two-dimensional harmonic oscillator have the 
explicit form: 



X + iSn-kV 



i-k\ 



s„ = sign(n) and L"{x) arc associated Lagucrre polyno- 
mials. 
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Consider now the hole states. A hole wave function 
is a complex conjugated electron wave function, and the 
hole charge is +e. Thus, we can obtain Hamiltonian of 
the hole in magnetic field from the electron Hamiltonian 
(2) by complex conjugation and reversal of the sign of 
the vector potential A(r2). In the representation of sub- 
lattices {^2,-62} it is 



Ho = 



H 



a2 

at 



(6) 



where the operators 02 = IhP2+ — *''2+/2Z_fr and = 
IhP2- + i''2-/2^// commute with oi, and obey the 
commutation relation [02, a^] = 1. Energies of the hole 
Landau levels are the same as these of electron Landau 
levels (3), but have an opposite sign. 

Hamiltonian of electron-hole pair without taking into 
account Landau level mixing is just the sum of (2) and 
(6), and can be represented in the combined basis of elec- 
tron and hole sublattices {A1A2, A1B2, B1A2, B1B2} as 



Ho — Hi + H2 



V2 







V 



02 Oi \ 

ai 

a2 

a+ a+ J 



(7) 



It is known [41] that for electron- hole pair in magnetic 
field there exists a conserving two-dimensional vector of 
magnetic momentum, equal in our gauge to 



Pi +P2 - 7r[H X (n -r2)] 

2c 



(8) 



and playing the role of a center-of-mass momentum. 
The magnetic momentum is a generator of simultane- 
ous translation in space and gauge transformation, pre- 
serving invariance of Hamiltonian of charged particles in 
magnetic field [42]. 

The magnetic momentum commutes with both 
the noninteracting Hamiltonian (7) and electron- hole 
Coulomb interaction V{ri — r2). Therefore, we can find 
a wave function of magnetoexciton as an eigenfunction 
of (8): 



*P«in2 (l"l,I"2) 



1 f „ [e^ X r] 

-exp|.R(p + 4^ 



x$„,„,(r- ro). (9) 



Here R = (ri -I- r2)/2, r = ri — r2, is a unit vector in 
the direction of the z axis. The wave function of relative 
motion of electron and hole $„j„2(r — Tq) is shifted on 
the vector ro = Ih[gz x P]. This shift can be attributed 
to electric field, appearing in the moving reference frame 
of magnetoexciton and pulling apart electron and hole. 

Transformation (9) from 4* to $ can be considered 
as a unitary transformation $ = U^, corresponding a 
to a switching from the laboratory reference frame to 
the magnetoexciton rest frame. Accordingly, we should 
transform operators as A — > SAS^. Transforming the 
operators in (7), we get: SaiS^ = &i, Sa^S^ = b'^ , 



Sa2S^ 



-62, Sb+S^ 



Here the operators 



bi 



IhP- - ir_/2lH, 



IhP+ + ir+/2lH, &2 = 
IhP+ ~ if+l'^^H^ ^2 = IhP- + ^f^/'^lH contain only 
the relative electron-hole coordinate and momentum and 
obey commutation relations [bijb'l] = 1, [&2,&2^] = 1 (all 
other commutators vanish). 

Thus, the Hamiltonian (7) of electron-hole pair in its 
center-of-mass reference frame takes the form 



^0 



Up 



V2 



I'M 



I 



\ 





K 




-b2 





bi 



-bi 



\ 

bi 

-b2 
/ 



(10) 



A four-component wave function of electron-hole relative 
motion <&nin2j being an eigenfunction of (10), can be con- 
structed by successive action of the raising operators b^ 
and 6^ (see also [20, 21]): 

l5ni,0+l5>i2,0— 2 



$„,„,(r) = {V2 
/s 



111 



V 



l,|n2|-l(l") 
Sni4'\ni\-l,\n2\i^) 
■Sn2'/'|ni|>2|-l(l") 



(11) 



<?^|«l||«2|(l") 



The bare energy of magnetoexciton in this state is a dif- 
ference between energies (3) of electron and hole Landau 
levels: 



(12) 



Here we label the state of relative motion by num- 
bers of Landau levels ni and n2 of electron and hole re- 
spectively. The whole wave function of magnetoexciton 
(9) is additionally labeled by the magnetic momentum 
P. In the case of integer filling, when all Landau lev- 
els up to v-th one are completely filled by electrons and 
all upper levels are empty, magnetoexciton states with 
ni > v, n2 < V are possible. For simplicity, we neglect 
Zeeman and valley splittings of electron states, leading 
to appearance of additional spin-flip and intervalley ex- 
citations [20, 21, 24]. 



II. INFLUENCE OF COULOMB INTERACTION 

Now we take into account the Coulomb interaction be- 
tween electron and hole V{r) = — e^/er, screened by sur- 
rounding dielectric medium with permittivity e. Upon 
switching into the electron-hole center-of-mass reference 
frame, it is transformed as V'{r) = V{r + ro). To ob- 
tain magnetoexciton energies with taking into account 
Coulomb interaction, we should find eigenvalues of the 
full Hamiltonian of relative motion H' = Hq + V' in 
the basis of the bare magnetoexcitonic states (11). As 
discussed in the Introduction, a relative strength of the 
Coulomb interaction is described by the dimensionless 
parameter 



e 



2.2 

e 



(13) 
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FIG. 1. Magnetoexciton dispersions _E„j„2(P), calculated in the first order in Coulomb interaction (dotted lines) and with 
taking into account mixing between 16 low- lying magnetoexciton states (solid lines). The dispersions are calculated at different 
filling factors v and different r^: (a) v = = 0.5, (b) v = = \, (c) = 0, Ts = 2, (d) v = \, = 0.5, (e) = 1, rs = 1, 
(f) V = \, = 2. Dispersions of 5 lowest-lying magnetoexciton states 712 — )■ ni, indicated near the corresponding curves, are 
shown. 



When e ^ 1, Tg ^ 1 and we can treat Coulomb inter- 
action as a weak perturbation and calculate magnetoex- 
citon energy in the first order in the interaction as: 

Due to spinor nature of electron wave functions in 
graphene, the correction (14) to the bare magnetoexci- 
ton energy (12) is a sum of four terms, each of them 
having a form of correction to magnetoexciton energy in 
usual two-dimensional electron gas [20-22]. Dependence 
of magnetoexciton energy on magnetic momentum P can 
be attributed to Coulomb interaction between electron 
and hole, separated by the average distance ro oc P . 
Calculations of magnetoexciton dispersions in the first 



order in Coulomb interaction (14) have been performed 
in several works [20-24]. However, such calculations are 
well-justified only at small enough rg, i.e. at large e. 
When £ ^ 1 (this is achievable in experiments with sus- 
pended graphene [43-40]), the role of virtual electron 
transitions between different Landau levels can be sig- 
nificant. 

To take into account Landau level mixing, we should 
perform diagonalization of full Hamiltonian of Coulomb 
interacting electrons in some basis of magnetoexcitonic 
states 'I'p„in2J where electron Landau levels n\> v are 
unoccupied and hole Landau levels <v are occupied. 
To obtain eigenvalues of the Hamiltonian, we need to 
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FIG. 2. Magnetoexciton energies at rest E^n^{P — 0), calculated with taking into account N electron and TV hole Landau 
levels, with stepwise increasing TV (crosses). The fits to these energies with inverse-square-root function (solid lines) and limiting 
values of E^^l^iP = 0) at TV — > cxd (dotted lines) are also shown. The results are presented for different filling factors i/ and 
different r^: (a) u ^ 0, = 0.5, (b) u = 0, Ts = I, (c) u ^ 0, = 2, (d) ly = 1, ^ 0.5, (e) = 1, rs = 1, (f) u ^ I, = 2. 



solve the equation: 



det 



0. 



(15) 



We can constrain our basis to TV^ terms, involving TV 
Landau levels for electron {ni = i/ + 1, . . . ,1/ + TV) and TV 
Landau levels for a hole (n2 = :^, — TV+1). Since the 
Hamiltonian commutes with magnetic momentum P, the 
procedure of diagonalization can be performed indepen- 
dently at different values of P, resulting in dispersions 
En^nziP) of magnetoexcitons, affected by a mixing be- 
tween TV electron and TV hole Landau levels. 

We present in Fig. 1 dispersion relations for 5 low- 
est magnetoexciton states, calculated with and without 



taking into account the mixing between 16 lowest-energy 
states. The results are shown for Landau level fillings 
ly = and ~ 1, and for different values of Vg. Close to 
P = 0, magnetoexciton can be described as a compos- 
ite particle with parabolic dispersion, characterized by 
some effective mass Af„j„2 = [(P En-i^n2iP) /dP'^]^^\p=o- 
At large P, the Coulomb interaction weakens and the dis- 
persions tend to the energies of one-particle excitations 
(12). However, the dispersion can be rather complicated 
structure with several minima and maxima at intermedi- 
ate momenta P ^ . 

We see that the mixing at small Tg has a weak effect on 
the dispersions (solid and dotted lines are very close in 
Fig. l(a.d)). However, at rg ^ 1 the mixing changes the 
dispersions significantly. We can observe avoided cross- 
ings between dispersions of different magnetoexcitons. 
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and even reversal of a sign of magnetoexciton effective 
masses (see Fig. l(b,c,e,f)). Also we see that the high 
levels are more strongly mixed than the low-lying ones. 
Similar results were presented in [20] for rg = 0.73 with 
conclusion that the mixing is weak. 

As we see, at large rg the mixing of several Landau lev- 
els already strongly changes magnetoexciton dispersions. 
Important question arises here: how many Landau lev- 
els should we take into account to achieve convergency of 
results? To answer this question, we perform diagonaliza- 
tion of the type (15), increasing step-by-step a quantity 
N of electron and hole Landau levels. For simplicity, we 
perform these calculations at P = only. Energies of 
magnetoexcitons at rest, renormalized by electron inter- 
actions due to breakdown of the Kohn theorem, are the 
most suitable to be observed in optical experiments. 

The results of such calculations of E^n.^{P ~ 0) as 
functions of N are shown in Fig. 2 by cross points. Wc 
found semi-analytically that eigenvalues of the Hamilto- 
nian under consideration should approach a dependence 



nin2 



nin2 



n 

^ n 



(16) 



at large iV. We fitted the numerical results by this de- 
pendence and thus were able to find the limiting values 

En-in-2 of magnetoexciton energies with infinite number number of solutions w 



using magnetoexciton wave functions (11) and energies 
(12) (see also [18, 32, 34-38]): 

n(g,a.) = 5 — ^^T^^"in.(g), (18) 
1 1 \ 

1 ^ , ,2 ^ 

(19) 
10 1/ 

where g = A is the degeneracy factor and /„ is the oc- 
cupation number for the ji-th Landau level, i.e. /„ = 1 
at n < v and /„ = at ji > (we neglect temperature 
effects since typical separation between Landau levels in 
graphcnc in quantizing magnetic field is of the order of 
room temperature [10]). The matrix between magne- 
toexcitonic wave functions in (19) ensures that electron 
and hole belong to the same sublattice, that is needed 
for Coulomb interaction in exchange channel treated as 
annihilation of electron and hole in one point of space 
and subsequent creation of electron-hole pair in another 
point. 

Unlike electron gas without magnetic field, having a 
single plasmon branch, Eqs. (17)-(19) give an infinite 



fir, 



of Landau levels taken into account. 

We see in Fig. 2 that the differences between mag- 
netoexciton energies calculated in the first order in 
Coulomb interaction (the crosses at iV = 1) and the ener- 
gies calculated with taking into account mixing between 
all Landau levels (dotted lines) arc very small at rg = 0.5 
(Fig. 2(a,b)), moderate at rg = 1 (Fig. 2(b,e)) and very 
large at rg = 2 (Fig. 2(c,f)). Since convergency of the 
inverse-square-root function is very slow, even the mix- 
ing of rather large (of the order of tens) number of Lan- 
dau levels is not suflScient to obtain reliable results for 
magnetoexciton energies, as clearly seen in the Fig. 2. 

Note that the mixing increases magnetoexciton bind- 
ing energies, similarly to results on magnetoexcitons in 
semiconductor quantum wells [47, 48]. 



III. MAGNETOPLASMONS 

Magnctoplasmons are collective excitations of electron 
gas in magnetic field, occurring as poles of density-to- 
density response function. In the random phase approx- 
imation, dispersion of magnctoplasmon is determined as 
a root of the equation 



J (g) , each of them can 
be attributed to specific inter-Landau level transition 
ri2 n-i, affected by Coulomb interaction [18, 37, 38]. 
Note that at g — >■ 0, when Coulomb interaction V{q) be- 
comes weak, dispersion of each magnctoplasmon branch 
712(9) tends to the corresponding single-particle exci- 
tation energy i!^n?ri2 • 

At Tg <C 1, wc can suppose that magnctoplasmon en- 
ergy 51„i„2 (9) does not differ significantly from the single- 
particle energy En-^n2 ■ In this case a dominant contribu- 
tion to the sum in (18) comes from the term with the 
given ni and ri2. Neglecting all other terms, we can write 
(18) as 



n(g,w) w g- 



E^°^ +iS' 



(20) 



and from (17) wc obtain an approximation to plasmon 
dispersion in the first order in the Coulomb interaction: 



^ni"2 (9) 



711712 



gV{q)Fny„^{q). 



(21) 



l-V{q)U{q,u;)^0, 



(17) 



where V{q) = 27re^/e(j is the two-dimensional Fourier 
transform of Coulomb interaction and Tl{q, w) is a po- 
larization operator (or polarizability). Polarization op- 
erator for graphcnc in magnetic field can be expressed 



Magnctoplasmons in graphene were considered with- 
out taking into account Landau level mixing in a man- 
ner of Eq. (21) in the works [20, 39]. Other authors 
[21, 24, 34] took into account several Landau levels, 
and the others [35-3N] performed full summation in the 
framework of the random phase approximation (17)-(19) 
to calculate magnctoplasmon dispersions. 

Here we state the question: how many Landau lev- 
els one should take into account to calculate magncto- 
plasmon spectrum with sufficient accuracy? To answer 
it, we performed calculations with successive taking into 
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FIG. 3. Magnetoplasmon energies ^In^n^t calculated in the lowest Landau level approximation (solid lines), with taking into 
account mixing between 2 (short dash lines) and 3 (long dash lines) Landau levels of electron and hole, and with taking into 
account mixing between all Landau levels (dotted lines). The results are presented for different filling factors i/ and different 
rsi (a) v = 0, Tb = 0.5, (b) u = 0, Tb = 1, (c) v — 0, rs — 2, (d) v = 1, rs = 0.5, (e) v = 1, rs = 1, (f) u — 1, rs — 2. Dispersions 
of 3 lowest-lying magnetoplasmon modes n2 ni, indicated near the corresponding curves, are shown. 



account increasing number of Landau levels at different 
v and Tg. In Fig. 3, dispersions of magnetoplasmons in 
graphene calculated numerically are shown. Results ob- 
tained without taking into account Landau level mixing, 
with taking into account a mixing of two or three lowest 
Landau levels and with taking into account all Landau 
levels are plotted with different line styles. 

As we see, even taking into account the mixing be- 
tween 2 Landau levels changes the dispersions consider- 
ably (see the differences between solid and short dash 
lines in Fig. 3). However, the calculations with mixing 
between 3 Landau levels (long dash lines) are already 
close to the exact results (dotted lines), except for the 



high-lying magnetoplasmon modes. It is also seen, that 
the mixing considerably changes the dispersions even at 
moderate Vg (see, e.g.. Fig. 3(d) at rg = 0.5). Note that 
the mixing usually decrease magnetoplasmon energies 
and does not affect the long-wavelength linear asymp- 
totics of their dispersions. 

Therefore we conclude here that convergency of mag- 
netoplasmon dispersions in rather fast upon increasing a 
number of Landau levels taken into account. Several low- 
est Landau levels are sufficient to obtain rather accurate 
results. On the other hand, calculations in the lowest 
Landau level approximation, i.e. without taking into ac- 
count the mixing, can give inaccurate results, especially 
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in a region of intermediate momenta q I 



IV. CONCLUSIONS 

We studied influence of Landau level mixing in 
graphene in quantizing magnetic field on properties of 
elementary excitations — magnetoexcitons and magne- 
toplasmons — in this system. Virtual transitions be- 
tween Landau levels, caused by Coulomb interaction, can 
change dispersions of the excitations in comparison with 
the lowest Landau level approximation. 

Strength of Coulomb interaction and thus a degree of 
Landau level mixing can be characterized by dimension- 
less parameter rg, dependent in the case of graphene only 
on dielectric permittivity of surrounding medium. By 
embedding graphene in different environments, one can 
change Ts from small values to rg « 2 [49]. 

We calculated dispersions of magnetoexcitons in 
graphene and showed that the mixing even between few 
Landau levels can change significantly the dispersion 
curves at rg > 1. However, at small rg the role of the 
mixing is negligible, in agreement with the other works 
[20, 25]. Then the question about convergency of such 
calculations upon increasing a number of involved Lan- 
dau levels have been raised. 

We performed calculations of magnetoexciton ener- 
gies at rest with taking into account stepwise increasing 
number of Landau levels and found their inverse-square- 
root asymptotics. By evaluating limiting values of these 
asymptotics, we calculated magnetoexciton energies with 
infinite number of Landau levels taken into account. We 
demonstrated that influence of remote Landau levels of 
magnetoexciton energies is strong, especially at large Tg. 
Also it was found that calculations with taking into ac- 
count even several Landau levels provide results, rather 
far from exact ones. 

Also dispersion relations of magnetoplasmons in 



graphene were calculated in the random phase approx- 
imation with taking into account different numbers of 
Landau levels. We showed that even few Landau lev- 
els for electron and hole are sufficient do obtain accurate 
results, however the lowest Landau level approximation 
(i.e. calculations without taking into account the mix- 
ing) provide inaccurate results, especially for intermedi- 
ate momenta and high-lying magnetoplasmon modes. 

In out study, we completely disregarded renormaliza- 
tion of single-particle energies due to exchange with filled 
Landau levels in the valence band of graphene, since this 
question was considered elsewhere [20, 21, 24, 30, 40]. In 
our article, we focused on the role of Coulomb interac- 
tion in the electron-hole channel only. In this context, 
an important result of our work is that breakdown of the 
Kohn theorem in graphene leads to strong corrections of 
magnetoexciton energies not only due to exchange self- 
energies, but also due to virtual transitions caused by 
Coulomb interaction between electron and hole. 

We considered magnetoexcitons in the ladder approx- 
imation and magnetoplasmons in the random phase ap- 
proximations without taking into account vertex correc- 
tions and screening. Estimating the role of these factors, 
especially in the strong- interacting regime at large Tg, is 
a difficult task and will be postponed for future studies. 

The results obtained in this work should be relevant for 
magneto-optical spectroscopy of graphene [28, 29, 31, 50- 
52] and for the problem of Bose-condensation of mag- 
netoexcitons [53-55]. Excitonic lines in optical absorp- 
tion or Raman spectra of graphene can give experimen- 
tal information about energies of elementary excitations. 
Magnetoexcitons and magnetoplasmons can be observed 
also as constituents of various hybrid modes — polaritons 
[56], trions [57], Bernstein modes [5n] or magnctophonon 
resonances [59]. 
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ident of Russian Federation for Young Scientists MK- 
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